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1. Introduction 



A linear polymer attached to a hard wall by an endpoint loses entropy due to the 
presence of the wall. The loss of entropy induces a repulsive force on the wall. These 
forces have been measured experimentally pQ [31 0] and decay with distance from the 
point where the polymer is attached. 

A simple model of the above is a directed path from the origin in the first 
quadrant of the square lattice. More precisely, let d+ be the number of directed 
paths from the origin given North-East and South-East steps in the half square lattice 
1? + = m) 6 1? | m > 0}. Since the path loses entropy due to the boundary of 
1? + (the X-axis) there is a net pressure on the X-axis (which we shall also call a hard 
wall). 

The pressure at a point x — {N, 0) in the X-axis may be computed by determining 
the number of paths which avoid the point x. Let this number be Then the 

loss of entropy is given by 

P n (a;)=logd+(a:)-logd+. (1) 

More generally, the paths may be interacting with the hard wall by adsorbing in it, or 
may be interacting in some other way. In that case there is a temperature dependence 
in the entropy loss, and the appropriate expression for the pressure on the hard wall 
is given by replacing d^(x) and d+ above by the corresponding free energies: 

P n (T; x) = T n (T; x) - F n (T). (2) 

The above approach has been used for a square lattice self-avoiding walk model 
of linear polymers grafted to a hard wall [5] . Exact enumeration of self-avoiding walks 
grafted to a hard wall showed that the pressure decreases with distance from the 
origin. In this paper we intend to add to these results by examining a directed version 
of this model, and to extend our results to include adsorbing walks, as well as a model 
of an adsorbing staircase polygon. 



1.1. The pressure due to directed paths 

The number of directed paths of n steps from the origin in the square lattice Z 2 is 
P n = 2". If the path is restricted to the half lattice Z?_, then it is a positive path, and 
if it is further constrained by having to end in the X-axis (which is the hard wall), 
then it is a Dyck path [5] . 

The number of positive paths of length n is given by 

(2n\ , {2n + 1\ 

T 2n = I n J and T 2n+1 = I ^ ^ J , (3) 

and the number of Dyck paths of length 2n is given by 

1 /2?i\ 

D 2n = —— ■ 4 

n + 1 \ n J 

From these expressions one may determine that the number of Dyck paths passing 
through the point x with coordinates x = (2N, 0). This is given by 

D 2N -D 2n _ 2N - {N + 1){n _ N + 1) [ N ) [ n _ N )■ (5) 



Thence, one may determine the pressure on the hard wall at the point x from equation 
(D to be 

PZ (2N) = log (D 2n - D 2N ■ D 2n -2N) - log D 2n = log (l - ° 2N '^- 2N ) . (6) 
Putting N = [an\ for some a €E (0, 1) and expanding asymptotically in n gives 

P 2 D n (2 LonJ) = — _ J__ + O (n" 5 / 2 ) . (7) 
wniv > ar'{l — a) J v 7 

A similar calculation may be done with positive paths. In this case the pressure 
on the point (2N, 0) is given by 

P 2 T „(2A0 = log (T 2n - D 2N ■ T 2n ^ 2N ) - logT 2 „ = log (l - D ™^-™ ^ . (8) 

Expanding asymptotically, and keeping the leading order term gives 

P? n (2an) = 1 + O (n- 5 ' 2 ) (9) 

■\/7rn 3 a 3 (l — a) V ' 

for pressure of positive paths of length 2n at a distance 2N = 2 [a n\ from the origin. 
1.2. Adsorbing paths and staircase polygons 

We shall generalise the results in equations (J7J and © to models of adsorbing Dyck 
paths, adsorbing directed paths, and a model of adsorbing staircase polygons grafted to 
the X-axis. While we extract the infinite n behaviour in each model, we will mostly 
be concerned with finite n behaviour, which we shall approximate by asymptotic 
expressions. 

A directed path with steps (1, —1) and (1, 1) from the origin in the half lattice 
7? + is illustrated in figure [Ha). Vertices of the path in the adsorbing line Y = are 
called visits, and they are weighted with the generating variable z (which is related to 
temperature by z = e -1 / T in lattice units). By convention, the origin, although itself 
a visit, is not weighted. A directed path with weighted visits is an adsorbing directed 
path. 

If the directed path is constrained to end in the adsorbing line, then it is a Dyck 
path (sec figure [2(b)). A Dyck path with weighted visits is an adsorbing Dyck path 

muni- 

The partition functions of adsorbing directed and Dyck paths are known to be 
given by 

L«/2J , v 

^)=E( rn/2 " +m )(-ir a") 

for adsorbing directed paths, and 

Dn[z) = y 4W+2 » ) {Z - I)™ (11) 

y ' ^ n + 2(m+ 1) V \n/2 \+mJ K ' y J 

m— ' 

for Dyck paths, see for example references [U[7]. 

The (extensive) free energies of these models are defined by 

Fl\z) =\ogT n {z), and (z) = logD B (z). (12) 
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a 



(b) 



Figure 1. Two models of an adsorbing polymer grafted to a hard wall, (a) A 
directed path from the origin in Z5j_. The path interacts with the X-axis with 
an activity z > 0. If z > 1 then the path is attracted to the X-axis, if z < 1 it 
is repelled. The path exerts a force on vertices of the X-axis, (b) This model is 
similar to (a), but both endpoints of the path are constrained to lie in the X-axis. 
This is a model of adsorbing Dyck paths. Observe that we use the convention 
that the vertex at the origin is not weighted by z. 



If (q, 0) is a vertex in the adsorbing line Y — 0, then define, similar to the 
above, the free energies of paths which avoids (q, 0), and denote these by T^{z; q) and 

The force (or pressure) on the vertex (q, 0) is given by the free energy differences 
P^(z;q)=^(z;q)- T T {z) and pr>( z . q )=j:D( z . q )-j:D( z ) (1 3 ) 

in the each model. Observe that Pn( z ; <z) = Pn{ z 'i (?) = if g is odd, so that a non-zero 
force can be obtained only for even values of q. 
For both models we show that for a £ (0, 1), 



lim P£(z;2 [an/2\) = lim F,f (z; 2 [an/2\) 



0. 



log(z - 1), 



if z < 2; 
if z > 2. 



(14) 



For finite values of n, we determine asymptotic approximations for P^(z;2 [an/2j) 
and Pn{z; 2 \_an/2\) (for n even). These are given by 



P^(z-2[an/2\) 



v /27m 5 o3(i - a) log 2 (z - 1) ' 
V2 



- log (z - 1) 



A T 



if z < 2 

if z = 2 
if z > 2 



12na(l-a)(z-l)(z-2) ! 
where 

A T = (a 2 - a + l)(z 4 + 8z 3 + 30z 2 - 32z + 16) - 6aV, 
for the model of directed paths. In the case of adsorbing Dyck paths, 

if z < 2 



(15) 



P°{z;2[an/2\) 



v/27m 3 a 3 (l - a) 3 log 2 (z - 1) ' 
V2 
ima(\ — a) 

- log (*-!)- 



< - 



.4 



D 



12na(l - a)(z - l)(z - 2)' 



if z = 2 
if z > 2 



(16) 
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for even values of n, where 

A D = (a 2 - a + l)(z 4 + 8z 3 + 30z 2 - 32z + 16). 

In the case of a model of adsorbing staircase polygons with both endpoints grafted 
to the adsorbing line (figure Efa)), a similar asymptotic analysis gives the following 
expressions for the pressure: 

5z 



P»{z;2[an/2\) 



2n(2-z)V™ 3 (i-a) 3 ' 
3 

2^Jima{l — a) 

3(4z 2 - 29z+ 16) 



if z < 2; 
if z = 2\ 



log (z - 1) 



tn(z - 2) 
3z 2 



8n(z-2) ; 



(17) 



if z > 2. 



I 2y/im 3 a 3 (l - a) 3 (z - l) 2 (z - 2) |log(z - 1)| 

2. The forces exerted by adsorbing directed paths in a half-space 

The partition function of adsorbing directed paths which passes through the vertex 
(q, 0) is given by D q (z)T n ^ q {z) (by equations (fTOjl and (|TTj) ). Thus the partition 
function of paths avoiding the vertex (<7,0) is T n {q) — D q (z)T n ^ q (z). Substitution 
into equation (|13p and simplifying gives 



P„ T (z; q) = log 1 



g,(z)T„_,,(z) 



(18) 



Observe that P^(z-q) is negative, since the pressure is directed onto the vertex (q,0) 
from above (and in the negative direction). 
A similar argument shows that 



P»(z;q) = log (l 



D q (z)D n - q (z) 



(19) 



in the case of adsorbing Dyck paths. 

A result in reference [2] shows that for z > the partition function D n (z) may 
be expressed in as 



D 2n (z) = 



where the C s = ^ (?) 
follows: C s < 4 s and thus 



z - 1 V z - 1 



1 OO 

-2) + -£C a 



1 



(20) 



are Catalan numbers. The summation can be bound as 



z - 1 



< 



4" 4 s 



z - 1 



4"z 2 



This shows that for z > 2 



< D 2n {z) < 



1 V z- 1 



z - 1 V z - 1 



(z-2) 2 

z(z - 1) ^4(z - 1) 
Jz~2y 



Since 4 < z 2 / (z — 1) if z > 2, the above proves that 



D n (z) 



(l + o(l)), ifz>2. 



(21) 



(22) 



(23) 



z - 1 VV^T, 

The result in equation (|2"3"| can be substituted into equation (flT)]) to give the 
following lemma: 



Lemma 1. The limiting force on the adsorbing line by a Dyck path if z > 2 is 
lim P°(z; 2 [an/2\) = - log(z - 1). 

n— >oo 

for any a G (0, 1), where the limit is taken through even values of n. □ 

A similar argument may be made for directed paths, using the following 
representation for the partition function 

T2n{z)=[^-A 9(z-2) + ^C a (l-a(z-2)) l^-j (24) 

^ ' s—n ^ ' 

from reference 2 . This shows that 



T n (z)=\^==j (l + o(l)), ifz>2. (25) 

Substitution of equations ([2^]) and ([23]) into equation (|T5|) shows that 
Lemma 2. TTie limiting force on the adsorbing line by a directed path if z > 2 is 
lim P^(z;2[an/2]) = - log(z - 1). 

n— »oo 

/or any 06 (0,1). □ 

Lemmas [1] and [5] show that the limiting force in the regime z > 2 is — log(z — 1) 
in both models. Observe that this force approaches zero as z — > 2 + . By 
equation (| 13[) the force cannot be positive, and thence lim„_ i . 00 Pn(z; 2 [a«/2j ) = 
linwoc Pj(z; 2 La^/2J) = if z G (0, 2]. 

Approximating the pressure 

The partition functions T n (z) and D n {z) will be approximated by using the Stirling 
approximation for the factorial 

n! = VW^e- (1 + ^ + O (£)) (26) 

to approximate the binomial coefficients. The summations will be approximated by 
an integral which itself will be estimated using a saddle point method. 

Substitute the Stirling approximation in the summands of equations (|10j) and 
take logarithms and simplify the results. Replace the summation over m by an 
integral. Fixing n, there is an m — [Sn\ which maximizes the summands, and where 
S G [0, 1/2]. The parameter S is determined by an asymptotic expansion in 1/n: Put 
n = 1/e, m = 8/e'm the summands, expand the result in e to leading order and 
simplify. This gives 

^(i^p^p'V^. (27) 

Take the derivative with respect to 6 and solve for <5 to find the saddle point at 

2 ol , (28) 



2z 

where we note that 5 cannot be negative. 

A similar approach in the case of T n (z) gives exactly the same result for the 
location of the saddle point, as expected. 



-0.25 




-0.75 



Figure 2. The pressure P„(z;q) for directed paths in for n = 128, z = 3/2 
and q = 2 [an/2\, plotted as a function of a. The magnitude of the pressure is 
large close to the origin and decays with increasing a towards the free endpoint 
of the path. The curve is the asymptotic expression for P^{z\q) as in equation 
(136 D . The dotted curve is the exact computed pressure, as determined from the 
partition functions in equations l |10|l and {TTJ . 



2.2. Case 1: 1 < z < 2. 

The general approach is to approximate the summations in equation (|10p and (jlip by 
integrals 



T n (z) 



[n/2\ 



/ y r™/2] + m 
for adsorbing directed paths, and 

;L«/2J 4m + 2 

D n (z) 



(z - l) m dm 



(z-l) m dm 



(29) 



(30) 



=0 n + 2(m + 1) V L"/ 2 J + m , 

for adsorbing Dyck paths (and where the binomial coefficients are generalised to 
continuous values of m by replacing factorials by Gamma functions). 

For asymptotic values of n, the integrals above are dominated by values of m close 
to the saddle points which are located close to m = by equation (|28[) and which have 
a spread of \Jn above m = 0. 

The approximation of T n (z) and D n (z) are achieved by substituting n — 1/e 2 
and m — a/e — a^/n in the summands above. Expand the results in e, and simplify. 

In the case of T„(z), the above prescription results in 



T n (z) 



(z-1) 



a /* e - 2a a da. 



(31) 



Evaluating the integral gives 

T n (z) ~ 2"- 1 e (" log2 ( z - 1 »/ 8 (l + erf (v 7 ^ \\og{z - . (32) 

An asymptotic expansion of this in n and keeping only leading order terms gives 

2 «+i 



T n {z) 



im |log(z — 1)| 



(33) 
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Figure 3. The pressure P®(z;q) for Dyck paths for n = 128, z = 3/2 and 
q = 2 \ an/2\, plotted as a function of a. The magnitude of the pressure is large 
close to the endpoints of the path, but is smaller for values of a towards the 
middle region of the path. The curve is the asymptotic expression for P^ ) (z; q) as 
in equation H37II . The dotted curve is the exact computed pressure, as determined 
from the partition functions in equation iJTTJ. 



The same procedure may be applied to D n (z): This shows that 



D n {z) 



(z-1) 



da. 



(34) 



Evaluating the integral produces 
D n (z) ~ 2 n+1 / 2 (-±= + |log(z - 1)| e "l°g 2 (*-i)/8 (i + erf (%/2^|log(z - 1)| /4 
Expanding this asymptotic in n, and simplifying gives 

2n+3 

D n (z) 



\/2im 3 log 2 (z - 1)' 
The above results produces the following approximations for the pressures: 

P„ T (z;2LW2j)^log( 1 



N /2 7 rn 3 a 3 (l - a) log 2 (z - 1) 



for directed paths, and 
P r f(z; 2LW2J ) ~ log [l- 



v /2vrn 3 a 3 (l - a) 3 log 2 (z - 1) 



(35) 



(36) 



(37) 



for Dyck paths. For large values of n the logarithms may be expanded, and the results 
should be compared to equations (O and ([5]) for z = 1. 

In figure [2] the approximation for (z;2 [an/2\) (solid curve) is compared with 
the exact values (dotted curve) for n = 128 and z = 3/2. Observe that asm oo, 
then P„(z; 2 \an/2\) approaches zero. 

Expanding the expression for the pressures asymptotically in n gives the 
expressions for z < 2 in equations (fT5|l and (fl6|) . 



-0.25 



-0.75 




Figure 4. The pressure P^(z;q) for directed paths in 7?, for n = 128, 



2 and 



q = 2 [an/2j, plotted as a function of a. The pressure is large close to the origin 
and decays with increasing a towards the free endpoint of the path. The curve 
and dotted points are determined from the exact expressions for the pressure, 
which is obtained from equations d 1 1) and Jill) on the one hand, and equation 
l|38ll on the other hand. 



2.3. Case 2: z = 2. 

The partition functions evaluate exactly in terms of factorials and Gamma functions 
if z = 2: 



T„(2) = T 



[(f) 

/Fr(*J2) 



i 



and D„(2) 



(n/2)!(n/2)! 



(38) 



Substituting these expressions into equations ([T8"|) and (fl"9"|) gives exact expressions 
for the pressure, which we do not reproduce here. Observe that as n — > oo, then it 
follows from the above that lim n _>oo P° ( 2 5 2 LW 2 J) = limn->«> P% ( 2 5 2 L«"-/ 2 J) = 0- 
The exact presssure Pj(2; 2 [an/2j) is plotted for n = 128 in figure ((4|). 

Asymptotic expressions for the pressure can be determined by using Stirling's 
approximation in equations (I38p . Substitution and simplification gives the follows 



P n T (2;2Lan/2j) = log 1 



P„ D (2;2Lan/2j)=log 1 



y27rna(i — a) 



2(1 - a - Vl^) 
Ttn-y/a (1 — a) 



+ O ( ri- 3 / 2 



(39) 
(40) 



for the pressure in these models. Expanding the pressures asymptotically in n gives 
the results for z = 2 in equations (fT5|) and (fT6|) . 



Case 3: z > 2. 



Asymptotic expressions in this regime are obtained by exploring the saddle point in the 
summands of D n (z) and T n (z). This saddle point is located at Sn in the asymptotic 
regime (where 5 is given by (z — 2)/2,z in equation ([28]) ) . 
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Figure 5. The pressure P„(z;q) for directed paths in for n = 256, z = 5/2 
and q = 2 [an/2\, plotted as a function of a. The pressure is large close to the 
origin and decays with increasing a towards the free cndpoint of the path. The 
curve is determined using the asymptotic expressions in equations H46H and (1476 . 
The dotted curve is the exact pressure (determined from equations IllOt and lilt ). 



The width of the saddle is proportional to y/n. Hence, use the Stirling 
approximation in the summands, put m = (^-^) n + cty/n and n = 1/e 2 , expand 
in e to O(l), and simplify. This gives the saddle point approximations 



1"°° „2 2 



. 2 2 



and 



Integrating a over (—00,00) gives the saddle point approximations 

These results are not unexpected, since substitution and simplification into equations 
(|18p and (|19|) give the results in lemmas Q] and O 

This also shows that the expansion to 0(1) in e does not produce expressions 
which give corrections for finite n effects. 

Finite n corrections to the pressures are obtained by improving the saddle point 
approximations above. Expanding to 0(e) instead gives the approximations 



where 



- 2 ~ 2 A 

-Atc 



_ (a 2 z 4 - (4a 2 - 3)^ 3 + (4a 2 + 9)z 2 - 18 z + 12) az 
At = 6(z-l) 2 (z-2) 
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Figure 6. The pressure P^ } (z;q) for Dyck paths in for n = 256, z = 5/2 
and q = 2 [an/2\, plotted as a function of a. The pressure is large close to the 
cndpoints of the path, and smaller towards the center. The curve is the asymptotic 
approximation of the pressure, obtained by using equation 1471 1. The dotted curve 
is the exact pressure (determined from equation ( 1 1011 and l llll l). 



and 

£\/n ( z — 2 

U n {z) ~ 

where 



a 2 z 3 - (2a 2 + 3)z 2 + 9 z - 6) az 



6(z-l) 2 

Expanding e"^ 6 = 1 - A T e + \A\<? + 0(e 3 ) and e- A ° e = 1 - A D e + \A 2 D e 2 + 0(e 3 ) 
in the above and integrating to a should give improved asymptotic expressions for 
T n (z) and D n (z) if z is not too large. Simplification gives 

and 

z 4 + 4(3n-2)z 3 -30(2n-l)z 2 + 16(2z-l)(3n-l) / z \" 
12n(z — l) z {z — 2) \y/z—l/ 

Substitution of the above in equations (|13|i and ()18|) and taking n — > oo again gives 
— log(z — 1), the expected limiting pressure. 

For finite values of n these results show a correction. In figure [5] a plot of 
P£(z,2 \ an/2\) against a is shown for n = 256 and z = 5/2. Similarly, in figure 
[6] a plot of the Dyck path case is shown by plotting P^(z, 2 \an/2\) as a function of 
a for n = 256 and for z = 5/2. 

3. The forces exerted by adsorbing staircase polygons 



A staircase polygon adsorbing in the upper half plane is illustrated in figure [71(a). If 
the left-most and right-most vertex in the staircase polygon is deleted, then it becomes 
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(a) 



(b) 



Figure 7. A model of staircase polygons, (a) A staircase polygon with its left- 
most and right-most endpoints grafted to the adsorbing line. If the first two and 
last two edges are deleted, then a pair of osculating directed paths are obtained. 
These paths avoid one another, and the bottom path interacts with the adsorbing 
line via the activity z. (b) A staircase polygon can be cut in a vertical line through 
one of the visits of the bottom path to the adsorbing line into a left and a right 
pair of paths. If the first two and last two edges of the polygon is removed, then 
each pair of paths consists of a Dyck path below a directed path. 



a pair of directed paths in the upper half plane, one path below the other, and avoiding 
each other. 

The bottom path is assumed to start in the origin, and by geometry the top 
path starts in the vertex with coordinates (0,2). In this situation we say that the 
staircase polygon is attached to the X-axis. Observe that if the final vertex in the 
bottom path has coordinates (X, Y), then the final vertex in the top path necessarily 
has coordinates (X, Y + 2). 

The bottom path may visit the X-axis, and these visits are weighted by z. By 
convention, the visit at the origin is not weighted. 

In this model we assume that the bottom path always ends in the X-axis in a 
point with coordinates (2n, 0) - these are grafted staircase polygons. Since we have 
reduced the model to two directed paths, we relax the conditions above and assume 
that the top path ends in a point with coordinates (2n,2j + 2). In this case the 
partition function of the model is known [2] . 

Hence, consider two directed paths in the square lattice, avoiding one another, 
and with steps (1,1) and (1,-1). Suppose the first path is a Dyck path starting at the 
origin and terminating in the vertex with coordinates (2n, 0). Suppose furthermore 
that visits of the first path to the X-axis are weighted by z. Let the second path start 
in the vertex with coordinates (0, 2) and terminate in the vertex with coordinates 
(2n, 2j + 2), as illustrated in figure [H Then the partition function of the pair of paths 
is 



mi— 7712—0 



where 

K(mi,m 2 ,j, n) 



(48) 



2?i + 3 W 2n + 3 \( z _ 1 )mi+^, 
n + m 2 + 2 J \ n + mi +j + 3 ' 



(2m 2 + l)(2mi + 2j + 3) (mi + m 2 + j + 2) (mi - m 2 + j + 1) 
{2n+ l)(2n + 2)(2n + 3) 2 ' 



13 





























































1 














< 


K 




ST 
-> 






















































/ 














s 


/ 










to 










\ 


































\ 






















< 












\ 


/ 






























































5 




5 


5 






5 


5 







Figure 8. A Dyck path below a direct path. By deleting the first two edges of 
this pair of paths, two directed paths which avoids one another are obtained. The 
bottom path is constrained to end in the adsorbing line; it is a Dyck path, and it 
adsorbs via the activity z in the adsorbing line. The top path is a directed path, 
and its endpoint has height 2j + 2 above the adsorbing line. 

One may check that C (0, z) = 1, Ci(0, z) = z, Ci(l, z) = z and C 2 (0, z) = z + 2z 2 . 
For a derivation of this result, see for example references [UES [7] but note the misprint 
in those expressions when compared to the above. Observe that the length of each 
path is 2n so that n is the half-length of each path. Similarly, the distance between 
the endpoints is 2j + 2, so that j + 1 is the half-distance between the endpoints. 

Putting j = gives the partition function of the model in figure EJa), minus the 
first two and last two edges. The expression above simplifies (see equation (4.188) in 
reference [2]) to 

k=0 v 7 

and the appearance of factors which are Catalan or ballot numbers should not be 
surprising. 

The pressure on a vertex with coordinates (q, 0) in the X-axis can be computed 
by determining the partition function of pairs of paths which avoids this vertex. This, 
in turn, can be done if one first determines the partition function of pairs of paths 
which passes through the vertex with coordinates (q, 0). This is illustrated in figure 
[7]^b), where the bottom path is constrained to pass through the marked visit. By 
cutting the polygon into two parts in the vertical line which passes through the visit, 
the pair of paths are divided into two sets of two paths each (one pair stepping from 
the left, and the other from the right), and each pair with endpoints a distance 2j + 2 
apart as illustrated in figure [8] (Reflect the paths on the right hand side of the cut to 
get it in the same orientation). 

In other words, the partition function of the model in figure [7](b) is given by 

n 

Z*(z;q)=J2C g (zd)Cn- q (zJ) (50) 

3=0 

The partition function of staircase polygons avoiding the visit (q, 0) is given by 

Z'(z;q)=C n (z,0)-Z*(z;q). (51) 



C n (0,z) = ^-( 2n ) 
n + 1 V n J 
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These give an expression for the net pressure on the vertex with coordinates (g,0): 



pS(z; q) = log(C„(z; 0)) - log(Z'(z; q)) = log 1 



(52) 



C n (z;0), 

In other words, to determine the force, one must compute both C n {z; 0) and Z*(z; q). 

In principle, P„ (z; q) can be determined exactly from the results above, but these 
expressions are complicated and not very informative. Hence, we shall approximate 
them to both estimate the pressure at finite values of n, and in the case that n — > oo. 
As in the case of directed walks, there will be different results for z > 2 (the adsorbed 
phase) and z < 2. 



3.1. z = 2: 

In the case that z = 2 and j = 0, equation (|4^)l simplifies to 



C„(2,0) 



2(2n + l) 



c- 



(53) 



where C„ = — tt f 2 ™) is Catalan's number. 

'* n+l V n / 

Summing j = to n instead in equation (j48|) gives 

^r(n + 2)- 



More generally, for fixed j, C n (2,j) simplifies to 



C n (2,j) 



2(2n+l)(j + l) 2 C„ 



2n 
n + j 



(54) 



(55) 



> + j + l)(n + j + 2)_ 

These results can be used to obtain exact expressions for the force P„ (2; g) in equation 
([52]). 

In particular, one may evaluate Z*(2; q) in equation (|50p exactly to 
2(n(3g + 2) - 3 9 2 + 1)I> + |)r(g + |)r(n - <? + ±) 16" 



f2" o) = 

V7r 3 r(g + 2)r(n + 3)r(n-g) 
Dividing by C„(2,0) and computing P^(2;q) in equation ([52]) gives 



P n S (2;g)=log 1- 



(n(3g + 2) - 3g 2 + 1) T(g + \)T{n -g + \)T 
2^T{q + 2)r(n - q + 2)T{n + |) 




(56) 



(57) 



If one put q — 2 [an/2\ wheren < a < 1 and take the limit n — > oo, then a symbolic 
computations program [6] shows that lim„_ i . 00 P„ (2; 2 [an/2j) = 0. This shows that 
the pressure is zero in the limit asn->oo and z = 2. 



3.2. Asymptotics 

By putting mi +m2 = k and summing over mi, equation (|48[) can be cast in the form 



Cn(Z,j)=2_^ 



k=0 



(n + l) 2 



2n + 2 
n + j + k + 3 



fc=0 



' (j + k + 2)(n - 2j(fc + 1) - 2fc - 1) 
_ 2(n + l) 2 (2n+ 1) 



(z-l) s 



2n + 2 



2n + 2 



n + fc + 2/ \n + j + 2 



(58) 
(z-l) k . 




_l K. , — — —I 

Figure 9. Exact pressures for adsorbing staircase polygons with z = 2, 
q = 2 [an/2] and for n £ {16,32,64,128,256}. The pressure is zero in the 
limit that n — ¥ oo. These points were computed by using equation J57t . which is 
the exact expression for P^(2\q). 



This form presents a single summand over k which one may consider for approximation. 

Isolate the summands of C n (z,j) above and cast them in terms of Gamma 
functions. This gives 

8(j + fc + 2)r 2 (n+|)(z-l) fc 16" 

1 n(n + l)(2n+l)T(n + j + k + 4)T(n~ j -k)' [ ' 

2(j + k + 2) (2j(k + 1) - n + 2k + 1) T 2 (2n + 2) (z - l) fc 

2 (2n+l)r(n + j + 3)r(n-j + l)r(n + fc + 3)r(n-/s + l)' 1 j 

Numerical evaluations of C\ and C2 for given n indicate that the summands are 
large for small k <C j (and j is also small) for z < 2, at both fc and j small for z = 2 
and at fe 3> j for z > 2 (and both k and j are large). In addition, both sums in 
the above make a substantial contribution of opposite signs, and so both must be 
examined to determine suitable asymptotics. 

Proceed by determining the dominant terms in the summands of equation (|58|) . 
The binomial coefficients will be approximated by using the Stirling approximation 
for factorials (see equation ([26]) ). Take logarithms of the first summand, substituting 
n = 1/e, j = S/e and k — a/e and expand the resulting expression in e. The leading 
term is 

'16(l-<5-a) 5+Q - 1 
(6 + a + iy 
Taking the derivative and solving for a gives 

a m = max{^-(5,0} (62) 

since a m cannot be negative. 

A similar treatment of the summand in the second sum in equation (|58jl gives the 
leading term in an expansion in e: 

'16(1- S) 5 - 1 ^ "^ a - 1 

(i + a )a+i( 1 + sy 
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Taking the derivative and solving for a gives 

a M =max{^,0} (64) 
since a m cannot be negative. 

3.2.1. z < 2: Proceed by approximating the summands C\ and C 2 above. Take 
logarithms of C\ and C 2l substitute the Stirling approximation for factorials (see 
equation ([26]) ) and simplify the results. Since the summands are dominated by k <C j, 
and j small, substitute j — S/e and n — 1/e 2 and expand to 0(e). Summing over k 
and replacing e = 1/ '^/n and S = j/y/n show that to leading order C\ and C 2 cancel. 

This shows that higher order terms must be determined in this case. Expanding 
to higher orders in e, summing over k and combining the contributions of C\ and C 2 
gives 

4(j + 1) (2j 2 (2 -z)+ 3(8 -z) + 6) e-^ 2 /" z 16" 
^(C 1+ C 2 )c ^5 (z _ 2) 4 ■ ( 65 ) 

One may extract the asymptotic behaviour for j — and j = 0(y/n). In the second 
case, the above can be simplified taking only the fastest growing terms in each factor. 
This shows that 

24^16" ., . n 

if j = 0; 



„ . 7rn 5 (2-z) 4 ' 

C n (z-j)=Y{S 1 + S 2 )^ { 1 . a/ (66) 
z — ' ' Szj e~ 3 ' 16™ 

^— rg— , ifj" = 0(vn). 

7rn & (2 — z) J 

The partition function in equation (150[) can be approximated from these last 
expressions, and similarly, one may approximate C n (z; 0) in equation (|49j) . This finally 
gives an approximation for P„ (z; q) in equation (|52p . 
In particular, one obtains that 

, f°° ( 8zf^nW 8zfe-?K n -ri 16(n-q) \ 

yC q (z;j)C n - q (z;3) ~ / g— -3 , -5— -3 dj. 67 

f~" Q Jo \ nq b {2-z) i n(n- q) J {2~ z) A J 

The integral can be readily done, and after division by the asymptotic expression for 
C„(z;0) one is left with the following approximation for the force: 

P^^-logf 1 j. (68) 

V 2^(q(n-q)f 2 {2-zy) 

Replacing q by 2 |_a 7T./2J and simplifying then gives 

F r f(z;2Lan/2j)~log(l 5z ). (69) 

7 Jy \ 2n^a 3 (l-a) 3 (2-z) 2 / V ' 

Hence, for a £ (0,1) P,f (2; 2 [a n/2j ) as n — > 00. 

In figure fTUl the pressure P^(z;2 \_an/2\) is plotted for a £ (0, 1), z — 3/2 and 
n = 256 comparing the approximate expression in equation (|69|) against the exact 
calculated pressure. 
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Figure 10. The pressure P„(z;q) for staircase polygons in 



for n 



256, 



z = | and q = 2 \ an/2\ , plotted as a function of a. The pressure is large close to 
the endpoints of the polygons. The curve is the approximate pressure (equation 
I|69[ l1 while the dotted curve are the exact values of the pressure. 



3.2.2. z = 2: The case z — 2 has been solved explicitly, as seen in equations 
(j5u) and (|57) . The expression for P^(2;q) in equation (j57) is in terms of Gamma 
functions, and these can be approximated when q = 2 [an/2\ for large n using 
Stirling's approximation. This shows that 



P„ s (2;2Lan/2j)=log 1 



2^/mra(l — a) 



o (V 3 / 2 ) 



(70) 



In figure [TT] a plot of the approximation and exactly calculated values of the pressure 
is presented. 



3.2.3. z > 2: The summands C\ and C-x in equations ([59]) and (|60|) provide the 
starting point. Numerical experimentation shows that C\ is dominated by terms with 
k = 0(n) and j = O (\fn) while j + k = 0(n) with a spread of the peak proportional 
to i/n. This is in particular confirmed by the result in equation (|62|) which indicates 
that the dominant values of j and k in equation (I59p are at k — ((z — 2)/ z)n — j and 

Hence, put n = 1/e 2 , k — (z — 2)/(ze 2 ) — 5/e + a/e and j = S/e in equation 
l|59p. Take logarithms and expand in e to O(l). Exponentiating and integrating the 
resulting expression gives the asymptotic expression 

(z - 2)4 n z 2n+1 H 9 

8n 



Ci(n,j) 



(1 



(1 



>(1))) 



(71) 



0F(z - l)»+J+3 n 3/2 

where the substitution S — j j \fn was made. Expanding to higher order in e 
before integrating gives subleading corrections, and does not alter this leading term 
correction. 

A similar approach should give an approximation to Ci. However, some care 
is needed in this case. Determining only leading term behaviour leads to incorrect 
results, and it is necessary to include higher order corrections. The arguments leading 
to equation (|64| show that the appropriate choices for j and k is n — 1/e 2 and 
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Figure 11. The pressure P^(z; q) for staircase polygons in for n = 256, z = 2 
and q = 2 \ an/2\, plotted as a function of a. The pressure is large close to the 
cndpoints of the polygons. The curve was determined from equation II57I I and the 
points along the dotted curve ares the exact values computed from the partition 
function in equation (58}. 



k = (z — 2)/(ze 2 ) + a/e and j = S/e. Substitute this into equation (fBT))) . and as 
before, take the logarithm, expand to 0(e) (rather than just 0(1)). Exponentiate the 
result and expand in e to 0(e 2 ). Integrating a € (—00,00) for small e then gives the 
asymptotic expression 



C 2 (n,j) 



z-2 

Jz~iy 



1 



z(2j + 1) - 4(j + 1) - 



2(^-2)j 2 



•o(l/n) 



(72) 



for C 2 . 

The approximation for C\ peaks sharply at j = while that for C2 dominates 
the contribution from C\ when j > and peaks at a value of j > 0. At j = both 
terms make a contribution. Hence, approximate the pressure at the point (g, 0) by 

■ (d (q, j) + C 2 (q, j)) (d (n - g, j) + C 2 (n- q, j)Y 



log 1 - 



(73) 



Ci(n,0) + C 2 (7i,0) 

Put g = 2 |_a 71/ 2 J and expand the integrand above and integrate term by term. 
Expand the results asymptotically in n and keep terms to 0{n~ 2 ). This finally gives 
the result 

1 3 4z 2 - 29z + 64 / 9z 



f£(z;2|an/2j)~log 



z - 1 



8 n(z - l)(z - 2) 
3z 2 



1 



8n(z-2) 2 



(74) 



2 A /7rn 3 a 3 (l-a) 3 (z - 2)(z - l) 3 log(z - 1) J 

Taking n — > 00 shows that P 7 f{z; 2 \ an/2\) — > — \og(z — 1), consistent with the result 
for adsorbing Dyck paths in the z > 2 regime (see lemma [2|). 
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Figure 12. The pressure P^{z;q) for staircase polygons in for n G 
{32,64,128,256}, z = 3 and q = 2 [an/2], plotted as a function of a. The 
pressure is large close to the endpoints of the polygons. The solid curve is the 
exact pressure at n = 256 (determined from equation 158t 1 while the set of dotted 
curves are the asymptotic results (equation 1741 ^ for n doubling from 32 to 256. 
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Figure 13. The pressure P^{z;q) for staircase polygons in for n G 
{32,64,128,256}, z = 4 and q = 2|an/2j, plotted as a function of a. The 
pressure is large close to the endpoints of the polygons. The solid curve is the 
exact pressure at n = 256 (determined from equation H58H *) while the set of dotted 
curves are the asymptotic results (equation 1741 ^ for n doubling from 32 to 256. 



4. Conclusions 



In this paper we have investigated the pressure exerted by a directed path in a 
half-space on the X-axes. This is a directed model of the forces exerted by a two 
dimensional polymer grafted to a hard wall. Our first model is of a directed path, in 
which we considered to two cases, namely an adsorbing Dyck path model with both 
endpoints grafted to the adsorbing line, and a directed path with only one endpoint 
grafted to the adsorbing line. 
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Figure 14. The pressure P^{z;q) for staircase polygons in for n G 
{32,64,128,256}, z = 5 and q = 2|an/2j, plotted as a function of a. The 
pressure is large close to the endpoints of the polygons. The solid curve is the 
exact pressure at n = 256 (determined from equation H58H *) while the set of dotted 
curves are the asymptotic results (equation 1741 ^ for n doubling from 32 to 256. 

The pressure curve for the Dyck path is both a function of the length of the path 
(n), and the adsorption activity (z). Asymptotic expressions were obtained for the 
pressure at (q, 0) = (2 [an/2j , 0) for a £ (0, 1). The pressure curve is symmetric about 
a = 1/2 and is asymptotically given by equation (|5T|) for z < 2, equation P0|) for z = 2 
and equation ([47)) for z > 2. The pressure is zero asn->oo and z < 2, but there is a 
net constant limiting pressure of magnitude log(z — 1) for z > 2. 

Similar observations can be made for the pressure due to an adsorbing directed 
path. In this model, asymptotic expressions for the pressure were determined and is 
given in equation (|37[) for z < 2, and by substituting equation (I39[) into equation (|19p 
for z = 2, and equation (|46p for z > 2. The pressure profiles for finite length paths 
are not symmetric in this model, but the limiting pressure is equal to the limiting 
pressure of Dyck paths. 

We have also determined the pressure due to adsorbing staircase polygons which 
were grafted to the adsorbing line at both ends. Determining the limiting pressure 
involved more careful analysis, but the results are given by expressions which are 
similar to the Dyck path results. The pressures are given asymptotically by equation 
([69)) for z < 2, equation dZQj) for z = 2 and equation ((74|) for z > 2. 
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